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Abstract 

It is shown that a smooth metal film (or a plasma layer) can be made transparent for 
an electromagnetic wave when two identical subwavelength diffraction gratings are 
placed on both sides of the film. The electromagnetic wave transmission through the 
metal film is caused by excitation of evanescent surface waves (plasmons) and their 
transformation into propagating waves at the gratings. A model which is developed 
analytically shows that the problem of the wave transmission is physically equivalent 
to the problem of excitation of two coupled resonators of evanescent waves which 
are formed at the two film surfaces. 
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1 Introduction 

The anomalous transmission of light through periodically perforated metal 
films has recently attracted enhanced attention due to the fundamental signif- 
icance of this effect and to its numerous possible applications, e. g. in optoelec- 
tronics, subwavelength lithography etc. It has been shown for the first time in 
Ref. 1 that a periodic set of subwavelength holes enhances at some frequencies 
the transparency of an optically thick film by several orders of magnitude. It 
was pointed out in this paper that this phenomenon can be related to reso- 
nance excitation of surface waves - plasmons - at the metal surface. These 
plasmons are trapped electromagnetic surface modes propagating along the 
interface between two media which have permittivities or permeabilities with 
opposite signs. The electromagnetic fields of the plasmons decay exponentially 
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into both media. Plasmons were discovered for the first time at the interface 
between a vacuum and a plasma-like medium [2] . 

Plasmons can be excited at a smooth dielectric-metal boundary by an incident 
electromagnetic wave when the wave is totally internally reflected at a near- 
boundary inhomogeneity of the dielectric permittivity [3,4,5,6]. An artificial 
periodical (e. g. corrugated or perforated) or natural accidental inhomogeneity 
of a metal boundary also allows the excitation of plasmons by an incident elec- 
tromagnetic wave. Wood's paper [7] of 1902 can be considered as the starting 
moment for investigations in this direction. In spite of such a long history this 
problem continues to attract attention (see, e. g., [8,9,10,11,12]). 

There have been many theoretical works which suggest models that try to clar- 
ify the physical mechanism that is responsible for the anomalous light trans- 
mission through the perforated metal film. Considerable part of these works is 
based on results of numerical simulations (see, e. g., [10,11,13,14,15,16,17,18,19,20,21,22]) 
and only some of them propose analytical investigations [9,24,25,26]. Most of 
the authors share the opinion that the enhanced transparency of the perfo- 
rated metal film is accompanied by excitation of plasmons. However, various 
authors consider different mechanisms of coupling between the plasmons on 
both surfaces of the film, which is the key element for the electromagnetic 
wave transmission. In Ref. [16], for example, the holes are considered as sub- 
wavelength cavities for evanescent waves coupling the plasmons on either side 
of the film. The theory developed in Ref. [15] considers coupling of the plas- 
mons through the evanescent electromagnetic modes in the holes which are 
treated as subwavelength waveguides. In Ref. [18] the enhanced transparency 
is explained by interference of diffracted evanescent waves. 

Other authors consider the holes in the film as a periodical inhomogeneity 
of the metal surface [10,19,20] or as a periodical modulation of the metal 
dielectric permittivity [24,9,25]. In these models plasmons on both sides of 
the film are coupled due to the overlapping of their fields in the metal body 
and periodical modulation of the surfaces or of the dielectric permittivity 
is necessary for coupling between incident and transmitted electromagnetic 
waves and plasmons. 

It is appropriate to mention here that there are also diametrically opposite 
points of view on the role of plasmons in the electromagnetic wave transmit- 
tance through the perforated (or corrugated) metal film. On the one hand, 
it is asserted in Ref. [21] that the excitation of plasmons decreases the trans- 
parency of the film. On the other hand, it is shown in Ref. [17] that excitation 
of surface waves, whatever its nature (spoof plasmons [23] or Brewster-Zennek 
modes [22], for example), and not only plasmons, leads to the increase of the 
film transparency. 
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We do not set ourselves to clarify the contribution of one or another physical 
mechanism in the transparency enhancement as an object of this work. In this 
paper we do present a simple model which allows us to investigate analyti- 
cally the plasmon mechanism of light transmission through an optically thick 
metal film. In order to exclude any other possible mechanism we consider a 
construction which is noticeably simple for analytical description and in which 
only the plasmon mechanism is possible. Instead of a perforated film, a smooth 
metal film with two diffraction grids placed on both sides of the film are con- 
sidered. Besides the considerable simplification of the analytical description 
and the clear physical treatment of the electrodynamical properties of this 
"sandwich" , such a construction may possess some technological merits. 

The metal film is an open (not bounded by walls) resonator of surface waves. 
In the absence of dissipation, the resonator Q-factor is infinitely large. The 
resonator excitation by an external source results in a significant growth of an 
eigenmode amplitude. Since the eigenmode fields are concentrated near both 
sides of the film, the eigenmode amplitude growth leads to an increase of the 
fields at the opposite side of the film with respect to the source position that 
can be regarded as an anomalous film transparency. The aim of this work is 
to show that such interpretation of the anomalous transparency is not only a 
convenient model, but also reflects adequately the physics of the process. The 
resonator properties - eigenfrequencies and geometry of eigenmodes fields - 
completely, qualitatively and quantitatively, determine the dependence of the 
film transparency on the system parameters. 



2 How a metal film can be made transparent for light? 

Let us start from the simplest problem considering a plane monochromatic 
electromagnetic wave incidence at a vacuum-metal interface. The electromag- 
netic properties of the metal will be described by its dielectric permittivity 
(Drude model) e p = 1 — u p /u 2 , where ui p and u are the metal electron plasma 
frequency and the incident wave frequency, respectively. In this approxima- 
tion there is no difference between a metal film and a plasma layer and we 
will use both these terms as synonyms. Let us assume that the z-axis is di- 
rected along the interface normal, the x-axis lies in the plane of incidence 
and the metal occupies a right-hand half-space. The incident wave fields are 
proportional to exp(ik x x + ik z z — icut), where k is the wave vector. The re- 
flected and the refracted wave fields are proportional to exp(ik x x — ik z z — iujt) 
and exp(ik x x + ik pz z — iujt) correspondingly. Here k pz = ^ k 2 e p — k 2 is the z- 
component of the wave vector in the metal and k = cu/c. Dissipation processes 
in the metal will be neglected. 

When the wave frequency uj is smaller than the plasma frequency u < u p , 
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then the ^-component of the wave vector in the plasma, k pz , is imaginary for 
an arbitrary value of k x , i. e. the refracted wave decays exponentially deep 
into the plasma and its amplitude decreases as exp(—\k pz \z). If the plasma 
occupies a whole half-space z > 0, then the incident wave energy is completely 
transmitted to the reflected wave energy. If the plasma layer thickness h p , is 
a finite quantity, then the wave passes partially through the plasma but the 
transmitted wave amplitude is exponentially small, ~ exp(— \k pz \h p ). 

The plasma permittivity is negative, e p < 0, when uj < uj p . As is well known, 
along the interface between two media in which the dielectric permittivities 
have opposite signs, a surface electromagnetic wave can propagate. The surface 
wave fields decay exponentially on either side of the interface. The magnetic 
field of this wave is parallel to the interface plane. Let the y-axis be directed 
along the magnetic field of the wave. The dispersion relation for the surface 
wave propagating along the vacuum-plasma boundary can be written in the 
form: 

\Jkl\e p \ + k 2 x ^ ^ 
\ £ v\\J^x ~ 



It is not difficult to show that Eq. (1) possesses a real solution k x (u) if uj 1 < 
ujp/2. This solution corresponds to an undamped wave - surface plasmon - 
propagating along the interface. 

The surface wave phase velocity oo/k x is directed along the interface and is 
less than the speed of light c. As a result, this wave cannot be excited by an 
incident propagating electromagnetic wave. The surface wave can be excited by 
an "incident" nonpropagating (evanescent) wave for which k x > k\. In this case 
the term "incident" means the wave whose amplitude decays exponentially as 



~ exp(— \Jk x — k^z) on approaching the interface from a source placed in the 
half-space z < 0. 

Let us assume that the magnetic field of the evanescent incident wave is di- 
rected along the y-axis and is equal to Hi exp(— kz), and the magnetic fields 
of the surface wave on either side of the interface, in the z < and z > 
half-spaces, are equal to H 2 exp(+kz) and H 3 exp(— k p z), respectively. Here 



k = yk^ — kQ, k p = \Jk%+ \e p \kQ and a common factor exp(ik x x — iut) is 
omitted. Using the standard boundary conditions it is easy to express the 
amplitudes H 2 and H 3 in terms of Hi. 

l + kp/(k\e p \) 2 
1 l-k p /{k\e p \y 3 1 l-k p /{k\e p \y K> 



Equations (2) contain a resonance denominator whose zero coincides with a 
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root of the dispersion equation (1). The physical meaning of these formulae 
is very simple: the interface forms an original resonator whose eigenmode 
is the surface wave. In the model under consideration (dissipationless and 
unbounded in the x and y directions) the resonator energy loss is equal to 
zero. The incident wave Hi plays the role of an external force (a pump mode), 
therefore the amplitude of oscillations in the resonator increases without limit 
when the incident wave frequency approaches the resonator eigenfrequency. 

When plasma occupies the space between two planes z = and z = h p , 
then the second surface forms the same resonator. The resonators on either 
sides of the plasma layer are connected by fields which decay exponentially 
inward the plasma body. Therefore the resonators coupling coefficient is small 
(it is assumed that k p h p 3> 1). A coupling of two identical resonators leads to 
eigenfrequency splitting into two frequencies which are shifted proportionally 
to the coupling coefficient. These eigenfrequencies correspond to symmetric 
and antisymmetric eigenmodes relative to the layer center. 

Such a binary resonator can be excited in the same way as the single one, 
by an external nonpropagating electromagnetic wave. Solving the problem of 
such a wave passing through the plasma layer, one can obtain the expression 
for the magnetic field amplitude H tr of the transmitted wave on the other side 
of the layer, in the region z > h p : 

k\e p \ e- kh ? 

Htr ~ ~kp AT (3) 



where 



D = -A 1 - ^ ]V - \ (l + ^)V^. (4) 




The equation D = is the dispersion equation whose roots determine the 
eigenfrequencies of the coupled resonators. Thus, in this case too the wave 
field amplitude approaches infinity when the pump wave frequency reaches 
one of the two eigenfrequencies of the binary resonator. 

If it were not a conditional application of "incident wave" and "transmitted 
wave" conceptions to the evanescent waves, it would be possible to interpret 
expression (3) not only as wave transmission through the plasma layer, but 
also as wave "amplification" by the layer. The "perfect lens" effect [27] is 
accounted for exactly by such "amplification" of evanescent waves by a left- 
handed material slab [28,29,30,31,32]. 

In order to obtain real transmission of a propagating wave through the plasma 
layer it is necessary to transform (partially or completely) this wave into an 
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evanescent one and to make an inverse transformation behind the layer. Since, 
as it follows from Eq. (3), the evanescent wave amplitude can be very large, 
essential enhancement of the wave transmission is expected even if the trans- 
formation coefficient is small. 

A one- or two-dimensional diffraction grating is a simple "transformer" of 
propagating and evanescent waves one into another. A wave incident at the 
grating gives rise to a set of reflected and refracted waves whose wave vectors 
projections k± ntm at the grating plane differ from the projection k± of the 

— * 

wave vector of the incident wave by an integer number of grating vectors k gi2 : 

— * — » — » — » 

k± n ,m = k± + nk gi + mk g2 , n, m — 0, ±1, ±2, . . .. 

Secondary waves for which k^ n ,m < ho, are propagating waves, and waves 
for which A;_L„ im > k are evanescent waves. The last ones are localized near 
the grating plane and decay exponentially, moving away from the grating. 
Independently of the incident wave polarization there are, in the general case, 
such evanescent waves whose magnetic field have non-zero projection onto the 
grating plane. These waves can excite the surface wave resonator. 

Thus it is expected that the plasma layer transparency for an electromagnetic 
wave, whose frequency is smaller than the plasma frequency, can be enhanced 
essentially when on each side of the layer two diffraction gratings are placed. 
It will be shown below that the transparency of such a "sandwich" reaches 
100%. 



3 The diffraction grating model 

In the general case only one of the evanescent waves scattered on the grating 
can be used for the resonator excitation. The contribution of all other non- 
resonant waves in the interaction of the incident wave with the system under 
consideration is small and will be neglected below. Therefore, without loss 
of generality, we may consider the grating to be one-dimensional and regard 
the magnetic fields of the incident and scattered waves as lying in the grating 
plane and being perpendicular to the grating inverse wave vector. 

The following dimensionless variables will be used below: £ = xuj p /c, ( = 
zuj p /c, VL = u/ojp. The wave numbers (real or imaginary) of wave spatial 
harmonics will be normalized to u p /c and designated in the vacuum and in 
the plasma as q and k, respectively. 

Suppose that on each side of a plasma layer, which is situated between the 
planes ( = d g and ( = d g + d p , two diffraction gratings are placed in the 
planes ( = and ( = 2d g + d p (see Fig. 1). In order to avoid cumbersome 
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expressions we will assume that the gratings are identical and placed at the 
same distances d g from the plasma boundaries, and the grating wave vector 
q g is directed along the £ axis. 

Let us consider a plane electromagnetic wave passing through the grating- 
plasma-grating (GPG) system. It is assumed that the wave falls from left to 
right on the GPG system, the wave vector of the incident wave lies in the £, ( 
plane and the wave magnetic field is parallel to the gratings planes. 



o d 



d +d 2d+d„ 

g p s p 



Fig. 1. Plasma layer between two diffraction gratings. 



The problem of a transmission coefficient calculation for a system with one or 
another shape of grating lies out of the framework of this paper. Therefore we 
will use the following simple model of wave scattering at a diffraction grating. 
Let us consider a modified wave equation for the wave magnetic field: 

o2 (ft U 

-q^ + -q^ + ^ H + 5 (0 fao + 2/ii cos k 9 H = 0. (5) 



Eq. (5) describes the wave passing through a 5-shaped layer situated in the ( = 
plane. The layer permeability is varied harmonically in the ^-direction. Such 
a layer posesses all the properties of a diffraction grating and therefore will be 
used for the description of waves scattering at the grating. The harmonic rule 
of the permeability variation does not restrict the model generality, because, 
as it was mentioned above, the anomalous transparency of the plasma layer is 
connected with the excitation of a surface wave which is resonant with one of 
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the spatial harmonics of the wave passing through the grating. Therefore, only 
one harmonic of the permeability variation plays a key role. The contribution 
of all the other harmonics is negligibly small and will be neglected further. 

The parameters /x and \i\ describe different properties of the grating. The 
parameter \x§ determines the mean transparency of the grating. The parameter 
fii determines the coefficient referred above of the mutual transformation of 
the propagating and evanescent spatial harmonics of the waves. Since the 
last parameter plays the key role, we can suppose at first that /i = 0. In 
accordance with Flouqet theorem let us look for Eq. (5) solution in the form: 

#(C,0 = e^£#„(C)e m ^. (6) 



A boundary condition for the expansion coefficients H(() at the 5-layer follows 
from Eq. (6): 



dHn 



+0 



dHn 
dC 



H n (-0)=H n (+0), 
+ /H [H n - 1 (0) + H n+1 (0)]=0. 



(7) 



The conditions (7) and standard boundary conditions at the plasma surface 
are enough for the problem solution. 



4 Electromagnetic wave transmission through the GPG system: A 
wave theory 



Let the incident wave frequency be f2. The wave magnetic field has the form 
Hifie tq ^ +m ^ , where go — \J& 2 ~ <zf ■ F° r the sake of simplicity let us assume 
that the grating period is less than the wavelength, i. e. all the spatial harmon- 
ics with n 7^ in Eq. (6) are evanescent. Only two harmonics are of interest 
for the problem: the zero harmonic because the total energy flux behind the 
second grating is concentrated only in it, and one from the evanescent har- 
monics with n^O (let it will be n = +1 for definiteness) which can play the 
pump mode role for the surface waves resonator. The contribution of all the 
other harmonics will be neglected below. 

The system of waves and their harmonics which have to be taken into account 
for the transmission coefficient calculation, is shown in Fig. 2. The harmonic 
amplitudes will be marked by two indices. The first one corresponds to the 
wave number and the second one corresponds to the harmonic number. The 
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propagation direction of the corresponding wave is marked by an arrow. As the 
"propagation direction" for the evanescent wave, the direction of its amplitude 
decrease is chosen. The propagating waves are represented by horizontal lines 
and the evanescent waves are represented by conventional curves which fall 
along the direction of the wave amplitude decrease. 
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Fig. 2. System of waves and their harmonics. 

Using standard boundary conditions at the vacuum-plasma interface and the 
condition (7) at the gratings, let us exclude all internal waves and leave only 
equations which allow connecting the amplitudes of the reflected (#2,0 an d 
#2,1) and transmitted (Hg and #9,1) wave harmonics with the incident wave 
amplitude Hi . These equations are the following: 



#9,0 






+Sia 14: H 2tl 






«22#2,0 




+<5ia 2 4-#2,i 


= /^2#1,0 




5 2 a 32 H 2 fi 


+#9,1 


+«34#2,1 


= 5 2 p 3 Hifi 




+5 2 a 42 H 2fi 




+e 2qid °D H 2jl 


= S 2 p 4 H lj0 



The following notations are introduced here: 
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The function D is determined by Eq. (4). The coefficients and $ do not 
vanish at jii = 0. 

The matrix of coefficients in Eq. (8) has a block structure: two 2x2 blocks 
which do not contain the parameter fj,i are placed diagonally, and two anti- 
diagonal blocks which are proportional to Therefore for fii = the system 
(8) is split into two independent systems for propagating (with index 0) and 
evanescent (with index 1) harmonics. For — a non-trivial solution for the 
propagated harmonics is possible only in the presence of the external wave H lfi 
(coefficient a 2 2 7^ always). However a non-trivial solution for the evanescent 
waves exists even when the external wave is absent. This solution exists if the 
determinant of the lower 2x2 block is equal to zero, or: 



The condition (9) of course coincides with the dispersion equation (4). 

For small but finite \i\ <C 1 the system (8) determinant does not vanish at the 
roots of the dispersion equation (8) but becomes small, ~ n\. It means that at a 
frequency which is close to one of the surface waves resonator eigenfrequencies, 
the transmission coefficient defined as \H§<a/ H\<o\ can be anomalously large. 

Considering 5\ and 5 2 as small, one can easily derive from Eq. (8) the following 
expression for the amplitude iJ 90 of the wave transmitted through the GPG 
system: 



where A, P and R are certain combinations of the system (8) coefficients a^, 



5 Optically thick plasma layer 

The expressions for the Eq. (10) coefficients become essentially simpler when 
the plasma layer is thick, which is the case of the most interest. The coefficient 
A determines the amplitude of the wave transmitted through the plasma when 
the gratings are absent. This coefficient is exponentially small, A ~ e ~ K o d p <g; 
1. Here k = + 1 — f2 2 is the wave decay decrement in the plasma. When 
exp(— Kodp) 1, the expressions for the coefficients R and P can be written 



Do 



0. 



(9) 




(10) 



A- 
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as: 

R ~ -A^$-BlB 2 _e 2 ^e 2i P cos 2 /3, 

P ~ -A-^—B+B_e Kldp e ip cos 
qi\£ P \ 



where 



«1 



and 



/3 = g ^g + arctan — r^- 

901% 



Using the functions B±, expression (4) for D can be rewritten as 

D = B 2 _e Kldp - B 2 + e- Kldp . (11) 



For Hi <C 1 the amplitude H 9 is maximal when the wave frequency is close to 
one of the resonator eigenfrequencies Q res which are defined by the dispersion 
equation roots, D (Q res ) = 0. The condition exp(—K, d p ) <C 1 means that also 
exp(— Kidp) <C 1 and that the dispersion equation roots are exponentially close 
to the roots of the equation = 0. Therefore, one may put ?i|£ p |//Ci — 1 

everywhere, except for B_, and write the expression for the amplitude H s o in 
the form: 

ifJ 2 B 2_ eKl d Pe -2 qi d ge 2if3 cog 2 p 

'9,o — —r—„ : — ~ ^2— : — ; — —#i,o- (12) 



q qi \(B 2 - e - 2 ^ d v) - i^B^e^i^e^ cos 



Near the resonator eigenfrequency, which is now determined by the equation 
B 2 _ = exp(—2Kidp), Eq. (12) can be presented as follows: 



#9,0- 



B^{Vt res )e Kldp e' 2q ^e 2il3 cos 2 f3 



2-gg=- (tt 2 - ) _ iJ^ e -2 gi d ge ip C0S p 

dilf. es V res J qoqi r- 



±ie 



e -2 qi d ge 2if3 CQS 2p 



eff 



(VI 2 - n 2 es ) - i£ 2 //e - 2 «i^e^cos/3 



#1,0 = 
#1,0, 



(13) 
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where 



p 2 A*l 



2q oqi (dB_/dW res ) 



and the signs ± correspond to the two roots of the dispersion equation _B„ (fl re s) 
± exp(— Kid p ). 

It is easy to see that the module of expression (13) is maximal at the frequency 
which is determined by the condition Re (f2 2 — f^es) — iel^e 113 cos (3 = and 
reaches unity. In other words the transmission coefficient of an electromagnetic 
wave through the GPG system reaches 100% at frequencies which are close to 
the surface waves resonator eigenfrequencies, and it does not depend on the 
value of the coefficient fii. This coefficient determines only the bandwidth of 
the resonant transparency. 

The independence of the upper limit of the transmission coefficient from the 
value of the coefficient /Ji, if it is small enough, confirms the statement that a 
concrete shape of a grating does not play an essential part. Various gratings, 
one- or two-dimensional, differ one from another by the value of the coefficient 
fj,i, which describes a transformation of the incident propagating wave into the 
evanescent mode that excites the resonator. 

These results of the approximate expression (13) analysis are confirmed by 
a numerical solution of the system (8). The dependencies of the transmission 
coefficient K = \Hgfi/Hifi\ on the dimensionless frequency Q for various values 
of the parameter /ii are presented in Fig. 3. For convenience, the dependencies 
KiVt — Vt^s) are shown. Here are the roots of the dispersion equation (1) 
for the isolated resonator. As it follows from the dependencies presented in 
Fig. 3, the plasmons excitation can lead not only to the transmission coefficient 
enhancement but to its suppression too. For small (thick curve /ii = 0.02 
in Fig. 3) the transmission coefficient in the close vicinity of the resonance 
frequencies can be much smaller than the transmission coefficient of the iso- 
lated plasma layer. In the model being used the plasmons excitation is the 
only mechanism which is responsible for the layer transparency enhancement. 
Therefore the film transparency suppression in some frequency regions cannot 
be used as an argument against the plasmon mechanism of the transparency 
enhancement (see [21,18]). 

For small /ii the dependence K(Q) has the form of two reflection-symmetric 
Fano resonance profiles [33], that is caused by interference between the two 
terms in the right-hand side of Eq. 10. The first term describes a non-resonant 
transmission of the wave through the plasma layer, the second one describes 
resonance transmission caused by the plasmons excitation [34,35]. When 
increases, the Q-factor of the surface waves resonator decreases (see next sec- 
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Q-Q (0) 

res 

Fig. 3. Transmission coefficient K as a function of the normalized frequency Q. 
Layer thickness d p = 3.0. 

tion), which in turn leads to broadening of the resonance frequency bands. The 
overlapping of bandwidths of two resonances makes the role of non-resonance 
transmission negligibly small and the dependence K(Q) acquires the typical 
form of two coupled resonators. 



6 Electromagnetic wave transmission through the GPG system as 
a resonator excitation 

The fact that the peak value of the transmission coefficient for small \i\ <C 1 
does not depend on the value of this parameter, has a simple explanation. The 
Q-factor of the resonator of the surface waves is determined by the energy flux 
out of the resonator. This flux is completely concentrated in the propagating 
waves. These waves arise at the diffraction gratings due to transformation of 
the evanescent waves, which are the eigenmodes of the resonator, into the 
propagating waves. The transformation coefficient is proportional to The 
resonator field amplitude at the grating is proportional to ~ e ~^ d a . Thus, 
the energy flux out of the resonator is proportional to ~ n\e~ 2qida and the 
resonator Q-factor is inversely proportional to this value, Q^es ~ [i\e~ 2qid9 . 
Further, the resonator pumping is carried out by the non-propagating har- 
monic, which appears due to the external wave incidence at the first grating. 
This harmonic amplitude is proportional to ~ /iiH 10 . The efficiency of an 
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arbitrary resonator excitation by an external field ^ e xt(C) is defined by the 
projection {^ ext ^* e ^) of this field onto the resonator eigenmode field ^ res ((). 
In the case under consideration ty ext ~ /! 1 iJ lj0 e _ ' ?1< % ^ res (() ~ e qi ^ ds >), i. e. 

(V ext ** res ) ~ / pnH lfi e-^e^- d ^d( ~ ^H lfl e-^. 
o 



The steady-state amplitude A of the oscillations in a resonator with a finite 
Q-factor is determined by the expression: 

^0 



fi 2 - fi 2 es - % W r JQ r 



where Q and f2 res are the external field frequency and the resonator eigenfre- 
quency, respectively. In the case under consideration we obtain: 

A = H 71 ~ H 10 - 



where the coefficient a ~ 1. 

The field amplitude H 9 at the system output is determined by the product of 
the eigenmode amplitude at the grating surface, A e~ gid9 , and the coefficient 
/ii of the transformation of the evanescent wave into the propagating wave: 

fl »~ ff '«!l.-flL-»«rW.' (14) 



Eq. (14) has just the same structure as Eq. (13). The peak value of the trans- 
mission coefficient, as it follows from Eq. (14), is about unity and does not 
depend on the parameter fj,i. Note, that in the estimations presented above 
of the resonator Q-factor and the resonator pumping efficiency the plasma 
thickness does not appear. It means that nominally an arbitrary thick metal 
film (plasma layer) can be made transparent using the construction under 
consideration. 

The surface wave resonator properties are changed when diffraction gratings 
are placed near the plasma layer. Not only a finite Q-factor appears, as it 
was shown above, but the resonator eigenfrequencies are shifted as well. The 
frequency shift can be determined in the following way. 

The dispersion equation Eq. (4) follows from standard boundary conditions 
at a vacuum-plasma interface. At the same time, from two possible solutions 
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of the Maxwell equations in vacuum, only such a solution is chosen that de- 
creases exponentially upon leaving the plasma boundary. This condition can 
be written as a Sommerfeld's condition analogue for evanescent waves: 

dH/dC\ dg _ - Ql H\ dg _ = 0, 
dH /d(\ dg+dp+ o + ?i H \ dg+dp+ o = 0- (15) 



In the presence of the gratings the conditions (15) are no longer satisfied. The 
"medium" at the left and at the right from the plasma boundaries can now 
be characterized by the surface impedance Z and the boundary conditions 
analogous to Eq. (15) take the form: 

dH/d(\ dg _ o + iQZH\ dg _ o = 0, 
dH/d(\ dg+dp+0 - iQZ H\ dg+dp+0 = 0. (16) 



Let us calculate the impedance considering, for example, the left plasma 
boundary. As it is represented in Fig. (4), the wave (4, 1) which "leaks" from 
the plasma is scattered at the grating and additional waves appear. Neglecting 
the wave (3, 0) transmission through the plasma layer, it is possible to consider 
that this wave is almost reflected from the boundary and is transformed into 
the wave (4,0). Using Eq. (7), one can connect the amplitudes of the waves 
(4, 1) and (3, 1) at the plasma boundary: 

2 

H 31 = --^(l + e ^+^K><h) e~ 2qid °H A u (17) 



where at is the phase shift under the wave (3, 0) reflection. 

At the plasma boundary the total magnetic field H\ of a spatial harmonic, 
which is proportional to ~ exp[i(g^ + q g )£\, is equal to H 31 + H 41 . The 
tangential component of the electric field of this harmonic, E^, is equal 
to —(i/Vl)qi(H Ay i — i?3,i). Thus, as it follows from the impedance definition 
Z = —E^i/Hi (the surface normal is directed oppositely to the £-axis) and 
from Eq. (17): 



Qi 



a 2 

l _ ^Jz2_ e md g +ta/2 e -2q ldg 



cos(qod g + a/2) 



(18) 



We put here ^ C 1. 

The impedance (18) is a complex quantity. The imaginary part of the impedance 
is responsible for the eigenfrequency shift and the real part describes the res- 
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Fig. 4. Resonator boundary impedance calculation. 

onator loss. Since the impedance imaginary parts of the "medium" with and 
without grating differ from one another by the value 



Im (SZ) ~ J^ e ~ 2qid ° sm(2q d Q + a), 



the resonator eigenfrequency shift is proportional to ~ n\ for small 

If the impedance of the "medium" surrounding the plasma layer is known, the 
resonator properties can be defined more precisely than it was done before. 

Using the boundary conditions at the vacuum-plasma interface and the con- 
ditions (16), which can be rewritten in the form: 



dH/d(\ dg _ + q* H\ dg _ = 0, 
dH / d (\d g +d P +o ~ 9* H \d g +d P +o = °> 

one can derive the following dispersion equation: 



(19) 



D = [B+( qi + g,) - - g,)] \ B -(qi " <?*) ~ B+( qi + q*)} e~ K ^ - 

[B-(qi + q*) ~ B+( gi - g,)] [B+(qi - q*) - + I*)} ^ = 0(20) 
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Here 

q* = -iZn = qi - t^-e-^e^ cos/3*, £ = q d g + a/2. 



Assuming that the plasma layer is optically thick, K\d p ^> 1, and /xf <d 1, it 
is possible to put B+ ~ 1 in Eq. (20) and, keeping terms of order of at most 
Eq. (20) can be presented as follows: 



D ~ Aq{e 



2 K\d v 



i ±?_e 

9o qi 



2q 1 dg e if3* 



COS /3* 



(21) 



The expression in the square brackets coincides exactly with the denominator 
in Eq. (12). 

The dependence of the resonator eigenfrequencies on the grating transparency, 
which is characterized by the parameter /zo in Eq. (5), can be obtained in 
an easier way putting Hi = 0. Replacing the second equation in (7) by the 
boundary condition 



dH n 



+o 



~d( 



+ /i # n (0) = 0, 



(22) 



it is easy to obtain the following dispersion equation: 



J£Le- qid °B + +(l-ir- - I ' 



Ho 
2<Zi 



2<Zi V 2gJ 



(23) 



When the plasma layer is rather thick, exp(nid p ) 3> 1, then the roots of the 
dispersion equation (23) are close to the roots of the first preexponential factor 
and Eq. (23) can be simplified: 



2<7i \ 2gJ 



0. 



(24) 



Let us restrict ourselves to the case when the gratings are situated close to the 
plasma surfaces 1 , when exp(qid g ) ~ 1. Using the definitions of the functions 

1 When the gratings are situated rather far from the plasma surfaces so that q\d g > 
1, then an extra pair of roots appears. These additional roots do not disappear for 
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B±, let us write Eq. (24) in the form: 

SM ( i _ M = i. (25) 



It is easy to see that the roots fl res of the dispersion equation (25) decrease with 
the growth of Ho and reach zero for /x = q% + q g . For Ho > q^ + q g Eq. (25) has 
no real roots. This dependence Q res (Ho), which has been obtained numerically 
from Eq. (23), is presented in Fig. 5. 




Fig. 5. Dependence of the resonator eigenfrequencies on the gratings transparency 
(parameter Ho)- 

The resonator considered above is a resonator of evanescent waves, therefore 
for Hi = its Q-factor is infinitely large. For Hi ^ the transformation of the 
evanescent waves into the propagating ones makes the resonator Q-factor finite 
and allows the excitation of the resonator by an external propagating wave. As 
in the case for Ho = 0, the GPG system transparency can reach 100% when the 
incident wave frequency is close to one of the resonator eigenfrequencies fl res - 

Ho — > oo and in some range of fio variation, four roots can exist simultaneously. 
The transmission coefficient K is extremal in the roots locations. In the considered 
frequency range and for q\d g > 1 one more mechanism exists, which is not connected 
with the surface waves resonator and also leads to the plasma layer transparency. 
The overall picture of the wave transmission is more complex in this case and will 
not be considered here. 



18 



The dependence f2 res (/io) described above is fully confirmed by the numerical 
solution of the system (8). For example the dependencies K(Q) for several 
values of /i marked by the vertical lines in Fig. 5, are presented in Fig. 6. 




Fig. 6. Dependence of the transmission coefficient K on the wave frequency for 
several values of hq: (a) jiQ = 0.1, (b) /^o = 0.4, (c) /io = 0.5. 



Fig. 6 also demonstrates that the eigenfrequencies of the plasma layers which 
are either isolated (/zo = 0) or surrounded by gratings can differ greatly 
and this difference cannot be treated as an argument against the plasmon 
mechanism of the transparency enhancement, as it was done, for example, in 
Ref. [16]. 

In closing this section let us note that the above-mentioned equivalence be- 
tween electromagnetic wave transmission through a plasma layer and excita- 
tion of two coupled identical resonators allows one to determine how one or 
another factors unaccounted for in this paper affect the result. Thus, for ex- 
ample, dissipative losses into the plasma layer can be taken into account in the 
following simple way. Making use of the analogy between a resonator and an 
oscillator, it is easy to see that the coupling between two oscillators (between 
two resonators, or, in other words, between the electromagnetic fields on either 
side of the layer) is destroyed when the generalized coupling coefficient c coup i 
is small as compared with the generalized friction coefficient Cf ric , 

Ccoupi Cfri c . (26) 



19 



In the case considered here, the coupling coefficient is determined by over- 
lapping of the evanescent fields of eigenmodes of two resonators , c coup i oc 
d p e~ Kldp . The friction coefficient is determined by dissipation in the layer, 
Cf ric oc dplme p . Thus, the resonant transmission through the layer can be 
observed when the following condition is satisfied: 

lme p e Kldp < 1. (27) 



Note that just this condition restricts the superresolution of the "perfect lens" 
[36]. 



7 Conclusion 



The merit of the suggested model is that all other possible mechanisms (waveg- 
uide connection between the film surfaces through small holes, action of holes 
as subwavelength cavities for the evanescent waves, interference of diffracted 
evanescent waves, etc) of the electromagnetic wave transmission through a 
perforated or corrugated metal film are absent a priori. Excitation of the sur- 
face waves resonator is the only mechanism in this model which is responsible 
for the electromagnetic wave transmission through a metal film. The knowl- 
edge of the resonator properties is sufficient for determining how the GPG 
system parameters affect the frequency and incidence angle for which the film 
transparency is maximal. 

The resonator mechanism can be considered as one of possible competitive 
mechanisms of light transmission through a perforated metal film. A period- 
ically inhomogeneous film surface can be considered as a diffraction grating 
situated directly at the film surface, and holes can be considered as an ad- 
ditional channel of connection between the resonators at the two opposite 
film surfaces. The increase of the coupling between the resonators decreases 
the destuctive influence of dissipation on the wave transmission through the 
film. The analysis of the dependence of the resonator eigenfrequencies on the 
system parameters shows that the eigenfrequencies of a smooth film and the 
grating-plasma-grating system can be significantly different. Therefore even a 
strong deviation of the frequency, which corresponds to the transparency peak, 
from the plasma layer (metal film) eigenfrequency, cannot be unambiguously 
treated as a sign of another, different from the resonator, mechanism of elec- 
tromagnetic wave transmission through the metal film. 
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